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Abstract 

We consider a supersymmetric (SUSY) Grand Unified Theory (GUT) based on the gauge 
group Gps = 5?7(4)c x SU(2)l x S77(2)r, which incorporates non-minimal chaotic infla- 
tion, driven by a quartic potential associated with the Higgs fields involved in the spontaneous 
breaking of Gps- The inflationary model relies on renormalizable superpotential terms and 
does not lead to overproduction of magnetic monopoles. It is largely independent of the one- 
loop radiative corrections and can become consistent with the current observational data on 
the inflationary observables, with the symmetry breaking scale of Gps assuming its SUSY 
value. Within our model, the strong CP and the \i problems of the minimal SUSY standard 
model can be resolved via a Peccei-Quinn symmetry. Moreover baryogenesis occurs via non- 
thermal leptogenesis realized by the out-of-equilibrium decay of the right-handed neutrinos, 
which are produced by the inflaton's decay. We consider two versions of such a scenario, 
assuming that the inflaton can decay to the lightest or to the next-to-lightest right-handed 
neutrino. Both scenaria can become compatible with the constraints arising from the baryon 
asymmetry of the universe, the gravitino limit on the reheating temperature and the upper 
bound on the light neutrino masses, provided that the gravitino is somehow heavy. In the 
second scenario, extra restrictions from the SU(4)c GUT symmetry on the heaviest Dirac 
neutrino mass and the data on the atmospheric neutrino oscillations can be also met. 
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1 Introduction 

Non-minimal Higgs inflation (non-MHI) [1-3] - see also Ref. [4] - is an inflationary model of chaotic 
type which arises in the presence of a non-minimal coupling between a Higgs-inflaton field and the 
Ricci scalar curvature, 1Z. It has been shown that non-MHI based on a quartic potential with a quadratic 
non-minimal coupling to gravity can be realized in both a non-supersymmetric [1] and a sypersymmet- 
ric (SUSY) framework [3], provided that the coupling of the inflaton to 1Z is strong enough. In most 
of the existing models, the inflaton is identified with the Higgs field(s) of the Standard Model (SM) 
or the next-to-MSSM {Minimal SUSY SM) [3] - for non-minimal inflation driven by an inflaton other 
than the Higgs field see Ref. [5-9] for non-SUSY models and Ref. [10-12] for SUSY ones. 

On the other hand, SUSY GUTs arise as natural extensions of Physics beyond the MSSM. Within 
their framework, a number of challenges - such as the p problem, the generation of the observed 
baryon asymmetry of the universe (BAU) and the existence of tiny but non-zero neutrino masses - 
which the MSSM fails to address can be beautifully arranged. The achievement of gauge coupling 
unification within the MSSM suggests that the breaking of the GUT gauge symmetry group down to 
the SM one, G S m = SU(3) C x SU(2) L x U(1) Y , occurs at a scale M GU T ^ 2 • 10 16 GeV through 
some Higgs superfields. Therefore, the latter arise naturally as candidates for driving non-MHI - for 
earlier attempts within non-SUSY SU (5) GUT see Ref. [13]. In a such situation the GUT gauge group 
is already spontaneously broken during non-MHI through the non-zero values acquired by the relevant 
Higgs fields. Consequently, non-MHI does not lead to the production of topological defects. Moreover, 
the potential of non-MHI possesses a non-zero classical inclination and so, the inflationary dynamics 
is largely independent of the radiative corrections. As a consequence, the vacuum expectation values 
(v.e.vs) which the Higgs fields acquire at the end of non-MHI can be exactly equal to the values required 
by the unification of the MSSM gauge couplings. Finally, the predicted inflationary observables are 
consistent with the fitting [14] of the seven-year data of the Wilkinson Microwave Anisotropy Probe 
Satellite (WMAP7), baryon-acoustic-oscillations (BAO) and Hubble constant (Hq) data. 

These features are to be contrasted with the widely adopted models of standard SUSY hybrid in- 
flation (HI) [15], where the spontaneous breaking of the GUT gauge symmetry takes place at the end 
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of HI and, thus, topological defects are copiously formed [16] if they are predicted by the symmetry 
breaking. This is because, the standard SUSY HI is typically driven by a singlet field whereas the 
Higgs fields are confined to zero where the GUT symmetry is unbroken. Avoidance of cosmologically 
disastrous topological defects can be obtained within smooth [16, 17] or shifted [18, 19] HI by using 
either non-renormalizable [16, 19] or renormalizable [17, 18] superpotential terms, which generate sta- 
ble inflationary trajectories with non-zero values for the Higgs fields. Some of the latter constructions, 
though, are much more complicated than the simplest original model [15]. In the cases of standard [15] 
and shifted [18, 19] HI, radiative corrections play an important role in creating the slope of the infla- 
tionary potential and the v.e.vs of the Higgs fields turn out to be mostly lower than the GUT SUSY 
symmetry scale, since the relevant mass scale is constrained by the normalization of the curvature per- 
turbation [14]. Finally, all types of HI suffer from the problem of an enhanced (scalar) spectral index, 
n s , which turns out to be, mostly, well above the current data [14]. For several proposals aiming to 
improve on the latter shortcoming of SUSY HI see Ref. [20-24]. 

In this paper we present a model of non-MHI, adopting a SUSY GUT model based on the Pati- 
Salam (PS) gauge group G PS = SU(4) C x SU(2) L x S77(2) R . Note that SUSY PS GUT models 
are motivated [25] from recent D-brane constructions and can also arise [26] from the standard weakly 
coupled heterotic string. Employing only renormalizable superpotential terms, we then show that the 
model naturally leads to non-MHI within SUGRA avoiding thereby the overproduction of unwanted 
monopoles. Also the inflationary observables turn out to lie within the current data [14]. Our model 
possesses a number of other interesting features too. The /i-problem of the MSSM can be solved 
[27] via a Peccei-Quinn (PQ) symmetry, which also solves the strong CP problem, and the proton is 
practically stable. Light neutrinos acquire masses by the seesaw mechanism [29] and the BAU can 
be generated through primordial non-thermal [30] leptogenesis. We single out two cases according to 
whether the inflaton decays to the lightest [31] or the next-to-lightest [32] right-handed (RH) neutrino. 
In both cases the constraints arising from the gravitino (G) limit [33-35] on the reheating temperature 
and the BAU can be met provided that the masses of G lie in the multi-TeV region. On the other 
hand, the second scenario gives us the opportunity to combine the calculation of BAU with the present 
neutrino data [36] and the prediction of Gps for the masses of the fermions of the third generation. 

The plan of this paper is as follows. We present the basic ingredients - particle content and struc- 
ture of the superpotential and the Kahler potential- of our model in Sec. 2. In Sec. 3 we describe 
the inflationary potential, derive the inflationary observables and confront them with observations. In 
Sec. 4 we outline the two scenaria of non-thermal leptogenesis, exhibit the relevant imposed con- 
straints and restrict the parameters of our model for each scenario. Our conclusions are summarized 
in Sec. 5. Details concerning the derivation of the mass spectrum of the theory during non-MHI are 
arranged in Appendix A whereas effects of instant preheating potentially important for some values 
of the parameters are discussed in Appendix B. Throughout the text, we use natural units for Planck's 
and Boltzmann's constants and the speed of light (h = c = &b = 1); the subscript of type , \ denotes 
derivation with respect to (w.r.t) the field \ (e-g-> ,xx = <5 2 /<9x 2 ); charge conjugation is denoted by a 
star and log [In] stands for logarithm with basis 10 [e]. 

2 The Pati-Salam SUSY GUT Model 

We focus on a SUSY PS GUT model described in detail in Ref. [19] - see also Ref. [37]. The rep- 
resentations and the transformation properties of the various superfields contained in the model under 
Gps, their decomposition under Gsm, as well as their extra global charges are presented in Table 1. 
Recall that, in the PS GUT models, the SM hypercharge Qy is identified as the linear combination 

Qy = Qt$ + Q(b-l)/2 where Q T -s is the SU(2)r charge generated by = diag (1, —1) /2 and the 
Q(B-L) is the SU(4) C charge generated by T^> = diag (1, 1, 1, -3) /2^/fj. Here T™ with m = 1, 2, 3 
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are the 3 generators of SU(2)r and Tq with a = 1, 15 are the 15 generators of SU (4) q with 
normalizations Tr (T^) = 5 ab /2 and Tr (Tj£T£) = 5 mk /2, where Tr denotes trace of a matrix. 

The ith generation (i = 1, 2, 3) left-handed (LH) quark [lepton] superfields, Uj a and c^ a - where 
a = 1, 2, 3 is a color index - [e^ and z/j] are accommodated in a superfield Fj. The LH antiquark 
[antilepton] superfields u? a and d? a [e? and tf] are arranged in another superfield F?. These can be 
represented as 



( 



qn Qi2 Qi3 h and F i 



) 



/ C \ 

its 
I? 



with 



%a = [d ia -u ia ) , h = (e, -i/jj , 9 ? a = | j and Zf = | ^ j . (2.1) 

The gauge symmetry Gps can be spontaneously broken down to Gsm through the v.e.vs which the 
superfields 



H c 



Q c H2 

l H ) 



and H c 



QH2 



QH3 



1°h) 



with 



acquire in the directions v c H and v c H . The model also contains a gauge singlet S, which triggers the 
breaking of Gps, as well as an 5£7(4)c 6-plet G, which splits into g a and g a under Gsm and gives [26] 
superheavy masses to d c Ha and d c Ha . In particular, G can be represented by an antisymmetric 4x4 
matrix as follows 

q _ f^abcSc 9i\ ^ q = fe a bc5c 9l ^ 



I -5a OJ t ~9i 







(2.3) 



Here G is the dual tensor of G, defined by G IJ = e ijkl Gkl which transforms under SU{4)q as 
UqGUq. Also e IJKL = sjjkl [^abcl is the well-known antisymmetric tensor acting on the SU (4)c 
[SU(3)c] indices with £1234 = 1 [£123 = !]• In the simplest realization of this model [19,26], the 
electroweak doublets H u and H^, which couple to the up and down quark superfields respectively, are 
exclusively contained in the bidoublet superfield IH, which can be written as 



(2.4) 



In addition to Gps, the model possesses two global U(l) symmetries, namely a PQ and an R 
symmetry, as well as a discrete Z™ p symmetry ('matter parity') under which F, F c change sign. The 
last symmetry forbids undesirable mixings of F and IH and/or F c and H c . The imposed £7(1) R 
symmetry, £7(1)^, guarantees the linearity of the superpotential w.r.t the singlet S. Although S does 
not play the role of the inflaton in our case - in contrast to the case of HI - we explicitly checked 
that the presence of a quadratic S 2 term would lead to the violation of the stability of the inflationary 
trajectory. Finally the £7(1) PQ symmetry, £7(1)pq, assists us to generate the ^-term of the MSSM. 
Although this goal could be easily achieved [38] by coupling S to IH 2 and using the fact that S, after 
gravity-mediated SUSY breaking, develops a v.e.v, we here prefer to follow Ref. [19,27] imposing a 
PQ symmetry on the superpotential and introducing a pair of gauge singlet superfields P and P. The 
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Table 1: The representations, the transformations under Gps, the decompositions under Gsm as well as the 
extra global charges of the superfields of our model. HereUc £ ST/(4)c, Ul £ S77(2)l, Ur £ SJ7(2)r and 
T, f and * stand for the transpose, the hermitian conjugate and the complex conjugate of a matrix respectively. 



PQ breaking occurs at an intermediate scale through the v.e.vs of P, P, and the /i-term is generated 
via a non-renormalizable coupling of P and IH. We do not adopt here the resolution to the /^-problem 
suggested in Ref. [38], since it introduces a renormalizable term which creates (i) a decay channel of 
the inflaton which leads to a high reheating temperature in conflict with the G constraint and (ii) a 
tachyonic instability in the H u — system during non-MHI - as occurring in Ref. [11]. Lifting both 
shortcomings requires an unnaturally small value for the relevant coupling constant in our scenario (of 
order 10~ 6 or so), which is certainly undesirable. Following Ref. [19], we introduce into the scheme 
quartic (non-renormalizable) superpotential couplings of H c to Ff, which generate intermediate-scale 
masses for the vf and, thus, masses for the light neutrinos, i/j, via the seesaw mechanism. Moreover, 
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these couplings allow for the decay of the inflaton into RH neutrinos, uf, leading to a reheating tem- 
perature consistent with the G constraint with more or less natural values of the parameters. As shown 
finally in Ref. [19], the proton turns out to be practically stable in this model. 
The superpotential W of our model splits into three parts: 

W = Wmssm + ^pq + ^hps, (2-5) 

where Wmssm is the part of W which contains the usual terms - except for the /x term - of the MSSM, 
supplemented by Yukawa interactions among the left-handed leptons and v f: 



W MSSM = yijFiHF* = 
Vij (HdeUe) - H^eUv] + H d J eQ ia d c ja - H u T eQ ia u c ja ^ , with e = 



1 
-1 



(2.6) 



HereQj a = ^u ia <j ia j and Lj = ^ e^j are the i-th generation 5{/(2)l doublet LH quark and 
lepton superfields respectively. Summation over repeated color and generation indices is also assumed. 
Obviously the model predicts Yukawa unification at Mqut since the third family fermion masses 
originate from a unique term in the underlying GUT. It is shown [39,40] that exact Yukawa unification 
combined with non-universalities in the gaugino sector and/or the scalar sector can become consistent 
with a number of phenomenological and cosmological low-energy requirements. The present model 
can be augmented [37] with other Higgs fields so that H u and H d are not exclusively contained in IH 
but receive subdominant contributions from other representations too. As a consequence, a moderate 
violation of the exact Yukawa unification can be achieved, allowing for an acceptable low-energy 
phenomenology, even with universal boundary conditions for the soft SUSY breaking terms. However, 
here we prefer to work with the simplest version of the PS model. 

The second term in the right hand side (RHS) of Eq. (2.5), Wpq, is the part of W which is relevant 
for the spontaneous breaking of U (1)pq and the generation of the \i term of the MSSM. It is given by 

W ^ = A «>^r " x "m Jr ( meItfe ) ' (2 - 7) 

where Ms ~ 5 • 10 17 GeV is the String scale. The scalar potential, which is generated by the first 
term in the RHS of Eq. (2.7), after gravity-mediated SUSY breaking is studied in Ref. [19, 27]. For 
a suitable choice of parameters, the minimum lies at \ {P)\ = \(P)\ ~ sJm^^M^. Hence, the PQ 
symmetry breaking scale is of order y/my^Ms c±. (l0 10 — 10 11 ) GeV and the /i-term of the MSSM 
is generated from the second term of the RHS of Eq. (2.7) as follows: 

-A M ^Tr (meIH T e) = fiH d J eH u n ~ A M ^, (2.8) 

which is of the right magnitude if ~ (0.001 — 0.01). Let us note that Upq has an additional local 
minimum at P = P = 0, which is separated from the global PQ minimum by a sizable potential 
barrier, thus preventing transitions from the trivial to the PQ vacuum. Since this situation persists at all 
cosmic temperatures after reheating, we are obliged to assume that, after the termination of non-MHI, 
the system emerges with the appropriate combination of initial conditions so that it lies [41] in the PQ 
vacuum. 

Finally, the third term in the RHS of Eq. (2.5), VFrps, is the part of W which is relevant for 
non-MHI, the spontaneous breaking of Gps and the generation of intermediate Majorana [superheavy] 
masses for uf [d° H and d c H ]. It takes the form 

Wups = AS (H C H C - M| s ) + X H H cT GeH c + \ n H c GeH cl + A^ i^jS , (2.9) 
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where Mpg is a superheavy mass scale related to Mgut - see Sec. 3.2. The parameters A and Mps can 
be made positive by field redefinitions. It is worth emphasizing that our inflationary model is totally 
tied on renormalizable superpotential terms, contrary to the model of shifted HI [19], where a non- 
renormalizable term added in the RHS of Eq. (2.9) plays a crucial role in the inflationary dynamics. 

Suppressing henceforth the color indices, we can express Wrps m terms of the components of the 
various superfields. We find 

^hps = XS {y c H v c H + e c H e% + u c H u% + d c H d% - M 2 S ) 

- 2X H [y c H d c H - e c H u c H ) g c + 2\ H u c H d c H g c 

~ 2A 5 {u c H d c H - e c H u c H ) g c + 2\ R u c H d c H g c 

+ {e°H4 + d c H d1 - u c H uf - u c H uf) 2 /M s . (2.10) 

Let us note in passing that the combination of two [three] color-charged objects in a term involves a 
contraction of the color indices with the symmetric [antisymmetric] invariant tensor 5 a b [£abc]> e.g., 
u c H u c H = KbU c Ha u c Hh [u c H d c H g c = e &hc u c Ha d c Hh gf\. 

According to the general recipe [3, 10], the implementation of non-MHI within SUGRA requires 
the adoption of a frame function, Q, related to the Kahler potential, K, as follows 

n _ - 3e -«« . - 3 + ST + ™ + + ^ - - f*» + h . c0 , 

mp mp 2mp mp mp 2mp ' 

(2.11) 

where the complex scalar components of the superfields H C ,H C ,G and S are denoted by the same 
symbol and the coefficients ks and kn are taken real. It is clear from Eq. (2.11) that we adopt the 
standard quadratic non-minimal coupling for Higgs-inflaton, which respects the gauge and the global 
symmetries of the model. We also added the fifth term in the RHS of Eq. (2.11) in order to cure 
the tachyonic mass problem encountered in similar models [3, 10] - see Sec. 3.1. In terms of the 
components of the various fields, K in Eq. (2.11) reads 

K = -3m 2 In (l - ^ + k s ^ + ^ (v c H v c H + e c H e c H + u%u% + d%d% + h.c.)) 

(2.12) 

with cf) a = v%,VH,e c H ,e c H ,u c H ,u c H ,d c Hj d c H ,g c ,g c ,S and summation over the repeated Greek indices 
- e.g. a and (3 - is implied 

In the limit where mp tends to infinity, we can obtain the SUSY limit, Vsusy, of the SUGRA 
potential, Vhf - see Sec. 3. Assuming that the SM non-singlet components vanish, Vsusy turns out to 
be 

Vsusy = A 2 | (v c H v c H - M 2 S ) | 2 + \ 2 \S\ 2 {\u c H \ 2 + \v%\ 2 ) . (2.13) 

On the other hand, assuming minimal gauge kinetic functions, the D-term scalar potential Vhd of the 
PS Higgs fields takes the form 

2 15 2 2 ^ 2 

Vhd = Y £ (h c *T£H c1 - H^T^*H c ^j + 9 — (h c *T^H c1 - H c ^T^*H c ^j , (2.14) 

a=l m=l 

where g is the (unified) gauge coupling constant of Gps- In terms of the components of H c and H c , 
Vhd reads: 

[ ipH^n + ~ " C H*e C H ~ v c H e C £? ~ (" C H e C H ~ + « ~ ^H? 

+ (v%u c H + v c H u% - v%u c H - v%u%f - (y%u c H - v c H u% - v%u c H + v c H u%f 

(2.15) 



Vhd = 



,; ^| 2 -I^| 2 ) 2 + (I^| 2 -I^| 2 ) 2 + --- 



2 
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Here the first line includes contributions arising from the sum over the SU(2) R generators and T| 
in Eq. (2.14), which are the well-known Pauli matrices. The second line includes contributions from 
the sum over the S77(4) c generators T 7+2a [T 8+2a ] (a = 1,2,3) with 1/2 [1/2] and -i/2 [i/2] in 
the a4 [4a] entries respectively and zero everywhere else, and the third for T^ 5 and Tp\ The ellipsis 
represents terms including exclusively the SM non-singlet directions of H c and H c . Vanishing of 
the D-terms is achieved for \i) c H \ = \v c H \ with the other components of H c and H c frozen at zero. 
Restricting ourselves to the D-flat direction, from Vsusy in Eq. (2.13) we find that the SUSY vacuum 
lies at 

(S) ~ and \(u c H )\ = \(u c H )\ = M PS . (2.16) 

Therefore, Wrps leads to a spontaneous breaking of Gps. The same superpotential, Whps, a l so gives 
rise to a stage of non-MHI as analyzed in Sec. 3. Indeed, along the D-flat direction = 3> 
Mps and S = 0, Vsusy tends to a quartic potential and so, IVhps can be employed in conjunction 
with K in Eq. (2.11) for the realization of non-MHI along the lines of Ref. [10]. 

It should be mentioned that soft SUSY breaking and instanton effects explicitly break U(\)r x 
U (1)pq to 1,2 x ^6- The latter symmetry is spontaneously broken by (P) and (P). This would lead to a 
domain wall problem if the PQ transition took place after non-MHI. However, as we already mentioned 
above, [7(1)pq is assumed already broken before or during non-MHI. The final unbroken symmetry 
of the model is Gsm x ^™ p - 



3 The Inflationary Scenario 

We below outline the salient features of our inflationary scenario (Sec. 3.1) and then, we present its 
predictions in Sec. 3.3, calculating a number of observable quantities introduced in Sec. 3.2. 

3.1 Structure of the Inflationary Action 

Following the conventions of Ref. [11], we write the action of our model in the Jordan frame (JF) as 
follows: 

Sm = f d 4 x^ QmpQft + mp^S^d^ - V m + ■ ■ , (3.1) 

where Q a p = Q, and <p a are identified below Eq. (2.12). The ellipsis represents terms arising from 
the covariant derivatives D fl H c and D^H C , and terms that include the on-shell auxiliary axial-vector 
field [3, 11], which turn out to be irrelevant for our analysis below. Here Vm = Q 2 Vhf/9 + Vhd, with 
Vhf being the Einstein frame (EF) F-term SUGRA scalar potential, which is obtained from PVhps in 
Eq. (2.10) - without the last term of the RHS - and K in Eq. (2.12) by applying [42] 

with K^K^ = 5? and F a = W HPS ,^ + K^W H ps/mp. 

If we parameterize the SM neutral components of H c and H c by 

v c H = he w cos^/^2 and v c H = he iS sm6 u /V2, (3.3) 

we can easily deduce from Eq. (2.15) that a D-flat direction occurs at 

= e = 0, 9 V = tt/4 and e c H = e c H = u c H = u c H = d c H = d c H = g c = g c = 0. (3.4) 

Along this direction, Vhd in Eq. (2.15) vanishes and so, Vm = Vhf takes the form 

f> 4 x2 ( x l ~ 4m PS ) 2 . , , n 1 2 1 k H 

VHio = m P h ps; with f = -- = ! + c n xl and c n = -- + ~ (3.5) 
lb/ z 6 D 4 
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Here mps = Mp§/rrip and Xh = h/m-p. From Eq. (3.5), we can verify that for S> 1, Shi 
in Eq. (3.1) takes a form suitable for the realization of non-MHI: the terms in the ellipsis vanish, 
and more importantly Vhi develops a plateau since mps <C 1 - see Sec. 3.2. The (almost) constant 
potential energy density Vmo and the corresponding Hubble parameter i^Hio along the trajectory in 
Eq. (3.4) are given by 



2ui 



16/ 2 



A 2 ™ 4 ~ V 1/2 

^ and H mo = __JIIQ_ 
164 



Amp 



V^mp 4y/3cn 



(3.6) 



We next proceed to check the stability of the trajectory in Eq. (3.4) w.r.t the fluctuations of the 
various fields. To this end, we expand them in real and imaginary parts as follows 



X 



X\ + 1X2 



V2 



X 



X\ + 1X2 

~7T 



where X = e c H , u c H , d c H , g c and x = e,u,d,g. 



(3.7) 



respectively. Notice that the field S can be rotated to the real axis via a suitable R transformation. 
Since along the trajectory in Eq. (3.4) the various fields, X and X, are confined to zero, the radial 
parametrization employed in Eq. (3.3) is not convenient here. Performing a Weyl transformation as 
described in detail in Ref. [1 1], we obtain 



Shi 



Mk 
P 



where 



K K 
R K 



^ c TZ x h' K 



f + R and V m = V HF + V HD /f 



(3.8) 



In deriving this result, we take into account that f :X <C f t h with % = S 1 , x±, X2, x\ and X2, and keep 
only terms up to quadratic order in the fluctuations % and their derivatives. To canonically normalize 
the fields v c H and u c H , we have to diagonalize the matrix Mk- This can be realized via a similarity 
transformation involving an orthogonal matrix Ur as follows: 



U K M K U K = 6\ag (/,/), where f = f + 6$ l x 2 h and U K 



1 

V2 



1 

-1 



(3.9) 



By inserting 1 = UkU k = U k IIk on the left and the right of Mk, we bring the second term of the 
parenthesis in the RHS of Eq. (3.8) into the form 



(3.10) 



along the trajectory in Eq. (3.4). Here 9± = (Q ± 0) /\/2. Consequently, we can introduce the EF 
canonically normalized fields, h,8 + ,9-,9 u and x, as follows - cf. Ref. [3, 10] 



dh = VJ 

dh f V / ' P 



1 , 6 44 n Jh6 + n 

h — ^ — , (7+ — — — , ti- 



ke- 



V2 



2f 



77 and X = W <3 ' U) 



Taking into account the approximate expressions for h - where the dot denotes derivation w.r.t the 
cosmic time, t - J and the slow-roll parameters e, rf, which are displayed in Sec. 3.2, we can verify 

that, during a stage of slow-roll non- MHI, 0+ ~ Jh9 + /V2 since Jh ~ V6m P , 0_ ~ M-/y/2f 

and 9 V ~ hOv/^ff since h/^JJ ~ mp/y^, and finally % ~ x/VT- F° r the latter, the quantity 
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f I f 3 / 2 , involved in relating x to X> turns out to be negligibly small, since f / f 3 ^ 2 = f t hh/f 3 ^ 2 = 
— X^e\rj\mp /2^/3c-ji. Therefore the action in Eq. (3.8) takes the form 

Shi = J d 4 x^ ^- + £ d^d^ - V m j , with V m = V HF + ^ (3.12) 

where stands for /i, 9 + , 9-,9 u , x±,X2, xi,X2 and S. 

Having defined the canonically normalized scalar fields, we can proceed in investigating the sta- 
bility of the inflationary trajectory of Eq. (3.4). To this end, we expand Vhi in Eq. (3.12) to quadratic 
order in the fluctuations around the direction of Eq. (3.4), as we describe in detail in Sec. A.l. In 
Table 2 we list the eigenvalues of the masses-squared matrices M 2 ^ = {d 2 Vm/dxadx^ with 

Xa = 9 + , 9-, 9 V , Xi ; 2, xi,2 and S involved in the expansion of Vhi- We arrange our findings into 
three groups: the SM singlet sector, S — v c H — v c H , the sector with the u c H , u c H and the e c H ,e c H fields 
which are related with the broken generators of Gps and the sector with the d c H ,d c H and the g c ,g c 
fields. 

As we observe from the relevant eigenvalues of the mass-squared matrices, no instability - as the 
one found in Ref. [11] - arises in the spectrum. In particular, it is evident that ks > 1 assists us to 
achieve m| > - in accordance with the results of Ref. [10]. Moreover, the D-term contributions to 
m| and m|_ - proportional to the gauge coupling constant g ~ 0.7 - ensure the positivity of these 

masses squared. Finally the masses that the scalars d\p acquire from the second and third term of the 
RHS of Eq. (2.9) lead to the positivity of m~ for Xh of order unity. We have also numerically verified 
that the masses of the various scalars remain greater than the Hubble parameter during the last 50 — 60 
e-foldings of non-MHI, and so any inflationary perturbations of the fields other than the inflaton are 
safely eliminated. 

In Table 2 we also present the masses squared of the gauge bosons, chiral fermions and gauginos 
of the model along the direction of Eq. (3.4). The mass spectrum is necessary in order to calculate the 
one-loop radiative corrections. Let us stress here that the non-vanishing values of v c H and v c H trigger 
the spontaneous breaking of Gps to Gsm- In particular we have the following pattern of symmetry 
breaking 

SU(A) C x SU(2) R -> SU(3)c x U(1) Y - 

Therefore, 9 of the 18 generators of S77(4)c x SU(2)r are broken, leading to 9 Goldstone bosons 
which are "eaten" by the 9 gauge bosons which become massive. As a consequence, 36 degrees of 
freedom (d.o.f) of the spectrum before the spontaneous breaking (18 d.o.f corresponding to 8 complex 
scalars, u c H ^u c H& ,e c H and e c H , and 2 real scalars, 9 and 9, and 18 d.o.f corresponding to 9 massless 
gauge bosons, Aq — A}$, A^, A 2 ^ and A- 1 -) of Gps are redistributed as follows: 9 d.o.f are associated 
with the real propagating scalars (9+,Xi- and X2+ with x = n a and e) whereas the residual 9 d.o.f 
combine together with the 18 d.o.f of the initially massless gauge bosons to make massive the following 
combinations of them A^,A^ and A 1 - - see Sec. A.2 of Appendix A. 

From Table 2 we can deduce that the numbers of bosonic and fermionic d.o.f in each sector are 
equal. Indeed in the S — v c H — v c u sector, we obtain 10 bosonic d.o.f and 10 fermionic d.o.f. Note that 
we consider S as a complex field and we take into account the 1 d.o.f of h which is not perturbed in 
the expansion of Eq. (A.l). Similarly in the u c H — u c H and e c H — e c H [d° H — d c H and g c — g c ] sector 
we obtain 32 [24] bosonic d.o.f and an equal number of fermionic d.o.f. Note also that the spectrum 
contains a massless fermion which must be present due to the spontaneous SUSY breaking caused by 
the tree-level potential energy density in Eq. (3.6). 

The 8 Goldstone bosons, associated with the modes x± + and X2- with x = u a and e, are not exactly 
massless since Vhi,/i / - contrary to the situation of Ref. [19] where the direction with non vanishing 
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Fields 


Masses Squared 


ElGENSTATES 


The S -v c h -v c h Sector 


2 real scalars 

1 complex scalar 

2 gauge bosons 


ml = mlxl (2A 2 (x? - 6) + 15# 2 /) /24/ 2 
m | = A 2 m|,x 2 (1 + Qc n ) /12J 2 / 3 ~ 4^ 2 j 
m | = A 2 m 2 x 2 (12 + x 2 /) (6%/ - 1) /6/ 2 / 
= bg 2 mpx 2 h /8f, to 2 = 


0+ 

s 

A\AW 


4 Weyl fermions 
1 Majorana fermion 


= X 2 mpx 2 i /2ff 2 
m\ = hg 2 mpx 2 h /8f 
m 2 =0 




The u c Ha - u c Ha (a = l, 2, 3) and e c H - e c H Sectors 


2(3 + 1) real scalars 
2(3 + 1) gauge bosons 


<- = ™ P x 2 (A 2 (4 " 3) + 3 5 2 /) /12/ 2 

2 2 

mi = mi 

e— w— 

m± = 9 2 mpx 2 h /Af 


^2+, 

ei_, e 2+ 

/ia± 4± 
' R 


4(3 + 1) Weyl fermions 


m \ = g 2 mlx 2 h /Af 


A^,C,^ 


The d c Ha - d c Ha and gl - g c a (a = 1, 2, 3) Sectors 


3 • 8 real scalars 


ml = m 2 x 2 (A 2 x 2 + 24\ 2 R f) /24/ 2 
ml = m 2 x 2 (A 2 x 2 + 24X 2 H f) /24/ 2 

™| + = « + 4 ^/) /4/ 2 
m|_ = m 2 x 2 (A 2 (x 2 - 3) + 12A 2 ,/) /12/ 2 


^a ^a 
Pi, 52 

ja ja 
M l+, u 2- 

ja ja 
«1_, «2+ 


3 • 4 Weyl fermions 


m | §d = X H m P X h/f 
m \. = ^H m Wh/f 


^± 
^+ 



Table 2: The mass spectrum of our model along the inflationary trajectory of Eq. (3.4). To avoid very lengthy 
formulas we neglect terms proportional to rrip S and we assume \h — Xfj for the derivation of the masses of the 
scalars in the superhelds d c H and d° H . The various eigenstates are defined in Sec. 3. 1 and Appendix A. 



(vjj) corresponds to a minimum of the potential. These masses turn out to be m X Q = Ampx^/2/. On 
the contrary, the angular parametrization in Eq. (3.3) assists us to isolate the massless mode 0_, in 
agreement with the analysis of Ref. [3]. In computing below the one-loop radiative corrections, V rc , in 
our model, we do not treat the residual Goldstone bosons as independent fields, since their associated 
d.o.f are "eaten" by the massive gauge bosons - for a different point of view, see Ref. [50]. After 
all, as we expect and verified numerically, V TC has no significant effect on the inflationary dynamics 
and predictions, since the slope of the inflationary path is generated at the classical level - see the 
expressions for e and rj below - and no significant running of the relevant parameters occurs - contrary 
to the SM or next-to-MSSM non-MHI. Employing the well-known Coleman-Weinberg formula [44], 
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we find 



V TC - Vsu c H u c H + V u c H u%e c H e% + V d%d%,g^-. ( 3 - 13 ) 



where the individual contributions, coming from the corresponding sectors of Table 2, are given by 

I ( m~ m/fi m- 



2 Tfi^-~~ 2 

+ 3mlln^-4m 4 ? ln -£f± - 4mi ln ^ | , (3.14a) 

x A 2 V>sv+ A 2 A 2 ' 

v « c H* c H*k*H = ^2 ( 2m i- ln + 6m ± ln a 2 " " 8m ± ln A 2 " ) ' (3 - Mb) 

3 / m 2 rni m\ m 2 ~ 

m 2 ~ ml- 
- 4m% In— ^-4m 4 r ln — ^ I . (3.14c) 

V>9d A 2 ^gd A 2 ' 

Here A is a renormalization mass scale. Based on the action of Eq. (3.12) with Vm — Vmo+V TC , we can 
proceed to the analysis of non-MHI in the EF, employing the standard slow -roll approximation [45,46]. 
It can be shown [47] that the results calculated this way are the same as if we had calculated them using 
the non-minimally coupled scalar field in the JF. 

3.2 The Inflationary Observables - Requirements 

Under the assumption that there is a conventional cosmological evolution (see below) after non-MHI, 
the model parameters can be restricted, imposing the following requirements: 

3.2.1 According to the inflationary paradigm, the horizon and flatness problems of the stan- 
dard Big Bag cosmology can be successfully resolved provided that the number of e-foldings, TV*, 
that the scale fc* = 0.002/Mpc suffers during non-MHI takes a certain value, which depends on the 
details of the cosmological model. The required JV* can be easily derived [23], consistently with our 
assumption of a conventional post-inflationary evolution. In particular, we assume that non-MHI is 
followed successively by the following three epochs: (i) the decaying-inflaton dominated era which 
lasts until the reheating temperature T r h, (N) a radiation dominated epoch, with initial temperature T r h, 
which terminates at the matter-radiation equality, (iii) the matter dominated era until today. Employing 
standard methods [8,23], we can easily derive the required iV* for our model, with the result: 

On the other hand, JV* can be calculated via the relation 



K = \ f H * Si = \ [ H * dh J 2 ^-, (3.16) 
m P V Bih m vJh f Vm th 



where h* [h*] is the value of h [h] when fc* crosses the inflationary horizon. Also /if [/if] is the value 
of h [h] at the end of non-MHI determined, in the slow-roll approximation, by the condition - see e.g. 
Ref. [45,46]: 

max{e(/if), |r/(/if)|} = 1, where 
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"4 / V m,h \ m|/VHi,fe\ 4m£ (1 + 4c^mji s ) 



^hi 7 2J2 \ y ffl i " 34^ 



(3.17a) 



, ~ 2 ^him m p (Vm,hh V m ,hJ,h\ 4m 2 (l + 4 C7 em| s ) 
and 77 = mp — ^ — = — £ — ^ ^ ^- ~ -75 — , (3.17b) 

are the slow-roll parameters. Here we employ Eq. (3.6) and the following approximate relations: 

J - >/6^, V m , h - (! + 4c^m 2 s ) and F m ,^ ^ m| (l + 4c^m 2 s ) . (3 18) 

The numerical computation reveals that non-MHI terminates due to the violation of the ^criterion at a 
value of h equal to hf, which is calculated to be 



e(h t ) = 1 =► hf = (4/3) 1/4 m P ^ (l + 4c^m| s ) /c n . (3.19) 
Given that hf <C /i*, we can write /i* as a function of iV* as follows 

^ " ^ (i + 4 l"J s K ^ = 2mp ^ (1 + 4c ^ s)/3c "- (3 " 20) 

3.2.2 The power spectrum P% of the curvature perturbations generated by h at the pivot scale 
k* is to be confronted with the WMAP7 data [14], i.e. 



P]f = -^—r 9mCK l V2 = -Xr J^4 - 4.93 • 10-. (3.21) 
n 2^/3 7rm| 1^(^)1 27rm2 \ 6e(/i„) 

Note that since the scalars listed in Table 2 are massive enough during non-MHI, the curvature per- 
turbations generated by h are solely responsible for P-%. Substituting Eqs. (3.17a) and (3.20) into the 
relation above, we obtain 



P V2 K W ^ AJV. (3 22) 

16^/27rm2 (l + 4crcm| s ) 12^/2ttc tc 



Combining the last equality with Eq. (3.21), we find that A is to be proportional to c-r., for almost 
constant iV*. Indeed, we obtain 

A ~ 8.4 • W- 4 ttcti/N* => c n ~ 20925A for iV* ~ 55. (3.23) 

3.2.3 The (scalar) spectral index n s , its running a s , and the scalar-to-tensor ratio r must be 
consistent with the fitting [14] of the WMAP7, BAO and H data, i.e., 

(a) n s = 0.968 ± 0.024, (b) - 0.062 < a s < 0.018 and (c) r < 0.24 (3.24) 

at 95% confidence level (c.L). The observable quantities above can be estimated through the relations: 

n s = 1 - 6e* + 2r% ~ 1 - 2/jV*, (3.25a) 
a s = I (4^ - (n s - l) 2 ) - 2£* ~ -2£* ~ -2/iV 2 (3.25b) 
and r = 16e* ~ 12/iV 2 , (3.25c) 

where £ = m p Vni,h V m hfTh/^m = m P V^V,h I J + 277?. The variables with subscript * are evaluated 
at h = h* and Eqs. (3.17a) and (3.17b) have been employed. 
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3.2.4 The scale Mpg can be determined by requiring that the v.e.vs of the Higgs fields take 
the values dictated by the unification of the gauge couplings within the MSSM. Since the highest mass 
scale of the model - see Table 2 - in the SUSY vacuum, Eq. (2.16), is 



m ±0 = y/hj2f o9 \{v c H )\ with f = f ((h)) = 1 + 4c n m 2 PS (3.26) 

(recall that mpg = Mpg/mp) we can identify it with the unification scale Mguts i-e. 

[5gMps ^ V2M GVT m P 

m± = V 2-/T = MGUT MPS = Wml-2c n M^ (3 " 27) 

The requirement 5g 2 mp > 2c^Mq UT sets an upper bound on c-jz, which however can be significantly 
lowered if we impose the requirement of Sec. 3.2.1 - see below. When c-jz ranges within its allowed 
region, we take M PS ~ (1.81 - 2.2) • 10 16 GeV. 

3.2.5 For the realization of non-MHI , we assume that c-ji takes relatively large values - see 
e.g. Eq. (3.8). This assumption may [2,48] jeopardize the validity of the classical approximation, 
on which the analysis of the inflationary behavior is based. To avoid this inconsistency - which is 
rather questionable [10, 48] though - we have to check the hierarchy between the ultraviolet cut-off, 
A = nip/cTi, of the effective theory and the inflationary scale, which is represented by Vmih*) 1 ^ 4 or > 
less restrictively, by the corresponding Hubble parameter, = Vm{h*) l l 2 / l \^3mp. In particular, the 
validity of the effective theory implies [48] 

(a) Vm(K) 1/4 < A or (b) H* < A for (c) c n > 1. (3.28) 



3.3 Numerical Results 

As can be easily seen from the relevant expressions above, the inflationary dynamics of our model 
depends on the following parameters: 

A, X H , X h , k s , cn and T rh . 

Recall that we determine Mpg via Eq. (3.27) with g = 0.7. Our results are essentially independent of 
Xh, X h and ks, provided that we choose some relatively large values for these so as m|_, m| and 
m| in Table 2 are positive for A < 1. We therefore set Xh = X H = 0.5 and ks = 1 throughout our 
calculation. Finally T r h can be calculated self-consistently in our model as a function of the inflaton 
mass, mi and the mass Mjpc of the RH neutrino into which inflaton decays, and the unified Yukawa 
coupling constant 2/33 - see Sec. 4.1. However the inflationary predictions depend very weakly on T r h, 
because T r h appears in Eq. (3.15) through the one third of its logarithm, and consequently its variation 
upon some orders of magnitude has a minor impact on the required value of A* , which remains almost 
constant and close to 55. 

In our numerical code, we use as input parameters h*, Mjpc and cn- For every chosen c-ji > 1 and 

Mjpc, we restrict A and h* so that the conditions Eq. (3.15) and (3.21) are satisfied. In our numerical 

1/2 

calculations, we use the complete formulas for the slow-roll parameters and in Eqs. (3.17a), 
(3.17b) and (3.21) and not the approximate relations listed in Sec. 3.2 for the sake of presentation. Our 
results are displayed in Fig. 1, where we draw the allowed values of c-jz (solid line), T r h (dashed line) 
and the inflaton mass at the SUSY vacuum - see Sec. 4.1 - m\ (dot-dashed line) [/if (solid line) and /i* 
(dashed line)] versus A (a) [(b)] for A% c = 10 11 GeV. Note that the decay of the inflaton into a RH 
neutrino with the mass above is kinematically permitted, for the depicted A's. The lower bound of the 
depicted lines comes from the saturation of the Eq. (3.28c). The constraint of Eq. (3.28b) is satisfied 
along the various curves whereas Eq. (3.28a) is valid only along the gray and light gray segments of 
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Figure 1: The allowed by Eqs. (3.15), (3.21), (3.28b) and (3.28c) values of c n (solid line), T rh - 
given by Eq. (4. 7) - (dashed line) and mi (dot-dashed line) [hi (solid line) and h* (dashed line)] versus 
A (a) [(b)] forks = 1, A H = A# = 0.5, M IP c = 10 11 GeV and y 33 = 0.5. The light gray and gray 
segments denote values of the various quantities satisfying Eq. (3.28a) too, whereas along the light 
gray segments we obtain h* > mp. 



these. Along the light gray segments, though, we obtain h* > mp. The latter regions of parameter 
space are not necessarily excluded [49], since the energy density of the inflaton remains sub-Planckian 
and so, corrections from quantum gravity can still be assumed to be small. As c-ji increases beyond 
5.6 • 10 3 , Ac-jimps becomes larger than 1, TV* derived by Eq. (3.20) starts decreasing and therefore, 
non-MHI fails to fulfil the relevant requirement. All in all, we obtain 

l<c^<5.6-10 3 and 5 • 1(T 5 < A < 0.25 for 53.9 < N, < 55. (3.29) 

From Fig. l-(a), we can verify our analytical estimation in Eq. (3.23) according to which A is 
proportional to c^. On the other hand, the variation of hf and /i* as a function of c-ji - drawn in 
Fig. l-(b) - is consistent with Eqs. (3.19) and (3.20). Note that the inclusion of the term 4c-jiMp S in the 
numerators of these relations is crucial in order to obtain a reliable result for A > 0.13 or c-ji > 3 • 10 3 

- cf. Ref. [11]. Letting A or c-ji vary within its allowed region in Eq. (3.29), we obtain 

0-964 < n s < 0.965, -6.5 < < -6.2 and 4.2 > > 3.5. (3.30) 

Clearly, the predicted a s and r lie within the allowed ranges given in Eq. (3.24b) and Eq. (3.24c) 
respectively, whereas n s turns out to be impressively close to its central observationally favored value 

- see Eq. (3.24a) and cf. Ref. [10]. 



4 Non-Thermal Leptogenesis 

In this section, we specify how the SUSY inflationary scenario makes a transition to the radiation 
dominated era, and give an explanation of the origin of the observed BAU consistently with the G 
constraint. The main features of the post-inflationary evolution of our model are described in Sec. 4.1. 
In Sec. 4.2 we describe the additional constraints that we impose on our setting, and finally we delineate 
the allowed parameter space of our cosmological model in Sec. 4.3. 
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ElGENSTATES 


Masses 


ElGENSTATES 


Masses 


The S 




- v c H Sector 


The u c Ha - u c Ha and e c H - e c H Sectors 


5h 




= V2XM PS /(J)f 


uf_, u| + ,ei_, e 2 + 




9 V 




y/5/2f gMp S 


' R 


gMp S /^ 


0+ 




V2XM PS /J f 




gMps/Vfo 


S 




V^XMps/yTo 


Ar,^,^ 


gMps/y/fo 






V5/2/o<7M PS 


The - and 5a - ^ Sectors 


AW 







91,92 


2X H Mp S /Jfo 






V2XM PS /^/To 


^a ^a 


2X R M PS /^/h 


A~ L ,Vv- 




y/5/2f gM PS 




2X H M PS /Vh 


ah,^ + 







^a ^a 


2X n Mp S /^fo 


The vf 


Sector 




2X H M PS /^fo 




M iP c 


= 2A^cM2 s /M s Wo 


r 9 ,± 


2A^M PS /V7o 



Table 3: The mass spectrum of our model at the SUSY vacuum of Eq. (2.16). We use the abbreviations 
(J) = J (h — 2Mps) , fo = 1 + 4c7jmp S and fo = fo + 24c|jmp s ~ Jq. The various eigenstates and symbols 
are specified in Sec. 3.1 and Appendix A. 

4.1 The General Set-up 

When non-MHI is over, the inflaton continues to roll down towards the SUSY vacuum, Eq. (2.16). 
There is a brief stage of tachyonic preheating [51] which does not lead to significant particle produc- 
tion [52]. Soon after, as discussed in the Appendix B, the inflaton settles into a phase of damped 
oscillations initially around zero - where Vhio has a maximum - and then around one of the minima 
of Vhio- Whenever the inflaton passes through zero, particle production may occur creating mostly 
superheavy bosons via the mechanism of instant preheating [53]. This process becomes more efficient 
as A decreases, and further numerical investigation is required in order to check the viability of our lep- 
togenesis scenario detailed below. For this reason, we restrict to A's larger than 0.001, which ensures 
a less frequent passage of the inflaton through zero, weakening thereby the effects from instant pre- 
heating and other parametric resonance effects - see Appendix B. Intuitively the reason is that larger 
A's require larger c^'s, see Eq. (3.23), diminishing therefore /if given by Eq. (3.20), which sets the 
amplitude of the very first oscillations. 

Nonetheless the standard perturbative approach to the inflaton decay provides a very efficient decay 
rate. This is to be contrasted with the SM (or next-to-MSSM) non-MHI, where the consideration of 
non-perturbative effects is imperative [52] in order to obtain successful reheating. Namely, at the SUSY 
vacuum u H and v H acquire the v.e.vs shown in Eq. (2.16) giving rise to the mass spectrum presented in 
Table 3. Note that the masses of the various scalars - contrary to the masses of the fermions and gauge 
bosons - are not derived from the corresponding formulas listed in Table 2 with the naive replacement 
Xh = 2mps, since terms proportional to mp s are neglected there. In Table 3 we also show the mass, 
mi, of the (canonically normalized) inflaton Sh = (h — 2Mpg) /Jq and the masses Mipc of the RH 
neutrinos, uf, which play a crucial role in our leptogenesis scenario - we assume the existence of 
a term similar to the second one inside In of Eq. (2.12) for uf too. From Fig. 1 we notice that mi 
increases with A - as in the case of HI, cf. Ref. [19] - only for A < 0.0013 or c-jz < 30. For larger A's 
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(J) = J(h = 2Mpg) ranges from 3 to 90 and so m\ is kept independent of A and almost constant at 
the level of 10 13 GeV. Indeed, if we express Sh as a function of Sh through the relation 



^~J where J = yj I + \m\j\ ((h)) = y/l + 24(^m PS (4.1) 

is obtained by expanding J given in Eq. (3.11) at leading order in Xh, we find 

\/2AM PS Amp irr 4 m P ,o r ,^ 1(T 4 ? „ „ 

mi ~ ^— — — ~ — ~ ~ 3 • 10 13 GeV for A > ~ 1.3 • 1(T 3 (4.2) 

foJo 2V3c n 4.2^3 4.2\/6mps 

where we make use of Eq. (3.23) - note that /q ~ 1. Apart from some fields in the vf sector, Sh is 
the lightest among the massive particles listed in Table 3 for A given in Eq. (3.29), since ( J) S> 1 and 
g, Xh and X H are taken larger than A. Therefore perturbative decays of 5h into these massive particles 
are kinematically forbidden. For the same reason narrow parametric resonance [51] effects are absent. 
Also Sh can not decay via renormalizable interaction terms to SM particles. 

The various decay channels of the inflaton are mainly determined by the Lagrangian part contain- 
ing two fermions - see Eq. (A. 15). The inflaton can decay into a pair of RH neutrinos (z/f) through the 
following lagrangian terms: 

Ch* = - Ai„c ¥p.^-(l- 12 Clz m 2 PS ) Shufuf + h.c. . (4.3) 

From Eq. (4.3) we deduce that the decay of Sh into Vf is induced by two lagrangian terms. The first 
one originates exclusively from the non-renormalizable term of Eq. (2.9) - as in the case of a similar 
model in Ref. [19]. The second term is a higher order decay channel due to the SUGRA lagrangian - 
cf. Ref. [54]. The interaction in Eq. (4.3) gives rise to the following decay width 



**> = l k m 4 l - iJ 3r w,,h = (1 - lWs) • {4A> 

where Mjpc is the Majorana mass of the RH neutrino into which the inflaton can decay - see Table 3. 
In addition, it was [54] recently recognized that within SUGRA the inflaton can decay to the MSSM 
particles spontaneously - i.e., even without direct superpotential couplings - via non renormalizable 
interaction terms. For a typical trilinear superpotential term of the form W y = yXYZ, we obtain the 
effective interactions described by the langrangian part 

Ci y = Gycn^^-Sh (x^ z + Yij x ^ z + Z^ Y ) + h.c. , (4.5) 

where y is a Yukawa coupling constant, tpx , ipY and tpz are the chiral fermions associated with the su- 
perfields X, Y and Z, and whose scalar components are denoted with the superfield symbol. For these 
scalars a term similar to the second one inside In of Eq. (2.12) is assumed so as to obtain canonically 
normalized scalars and spinors through relations similar to the last equalities in Eqs. (3.11) and (A. 14) 
respectively. Taking into account the terms of Eq. (2.6) and the fact that the adopted SUSY GUT pre- 
dicts Yukawa unification for the 3rd generation at Mps, we conclude that the interaction above gives 
rise to the following 3 -body decay width 

P Uc % 3 3y 3 2 3 3/mA 2 M PS /o 3/2 
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with 1/33 ~ (0.4 — 0.6) being the common Yukawa coupling constant of the third generation computed 
at the mi scale and summation is taken over color and weak hypercharge d.o.f, in conjunction with the 
assumption that m\ < 2M 3 „c which is valid for both inflaton-decay scenaria considered below. 

Since the decay width of the produced v\ is much larger than T\- see below - the reheating tem- 
perature, T r h, is exclusively determined by the inflaton decay and is given by [55] 



/ 72 \ 1/4 , 

' rh = V5tt2 — J with Ti = IV + ri y , (4.7) 



where counts the effective number of relativistic degrees of freedom at temperature T r h- For the 
MSSM spectrum, g* ~ 228.75. From Fig. 1 we remark that T r h does not exclusively increase with A, 
but rather follows the behavior of m\. For M\y C = 10 11 GeV, we find that T\ y dominates over T\ u c for 
A > 0.002 or - via Eq. (3.23) - c n > 50. For 0.0013 < A < 0.03, T rh remains almost constant since 
/q ~ 1 - see Eq. (4.6). For A > 0.03, /q ~ 1 + 12c-^.mp S starts to deviate from unity and so, T r h 
increases as shown in Fig. 1. 

If T r h < Mjpc, the out-of-equilibrium condition [31] for the implementation of leptogenesis is 
automatically satisfied. Subsequently z/f decay into H u and L* via the tree-level couplings derived 
from the second term in the RHS of Eq. (2.6). Interference between tree-level and one-loop dia- 
grams generates a lepton-number asymmetry ej_, [31], when CP conservation is violated. The resulting 
lepton-number asymmetry after reheating can be partially converted through sphaleron effects into 
baryon-number asymmetry. However, the required T r h must be compatible with constraints for the G 
abundance, Y~, at the onset of nucleosynthesis (BBN). In particular, the B yield can be computed as 

(a) y B = -0.35y L with (b) Y L = 2e4^^. (4.8) 

4 1 1 m\ 

The numerical factor in the RHS of Eq. (4.8a) comes from the sphaleron effects, whereas the numer- 
ical factor (5/4) in the RHS of Eq. (4.8b) is due to the slightly different calculation [55] of T rh - 
cf. Ref. [31]. In the major part of our allowed parameter space - see Sec. 4.3 - T\ ~ Y\ y and so the 
involved branching ratio of the produced uf is given by 

rv ^ r^c _ vr 2 (i - i2c n m 2 PS ) 2 (4 9) 



Ti Tiy 72c^y| 3 m^ s m 



For Mipc ~ (10 11 — 10 13 ) GeV the ratio above takes adequately large values so that is sizable. 
Therefore, the presence of more than one inflaton decay channels does not invalidate the non-thermal 
leptogenesis scenario. On the other hand, the G yield at the onset of BBN is estimated to be [35]: 



Y d ~c 5 T lh with c s = 1.9 ■ 10 _2 7GeV. (4.10) 

Let us note that non-thermal G production within SUGRA is [54] also possible. However, we here 
prefer to adopt the conservative approach based on the estimation of Yq via Eq. (4.10) since the latter 
G production depends on the mechanism of the SUSY breaking. 

Both Eqs. (4.8) and (4.10) calculate the correct values of the B and G abundances provided that 
no entropy production occurs for T < T r h, as we already assumed - see Sec. 3.2. Regarding the 
estimation of el, appearing in Eq. (4.8), we single out two cases below, according to whether the 
inflaton can decay into the lightest (Pf) or to the next- to-lightest (P|) RH neutrino. Note that the decay 
of the inflaton to the heaviest RH neutrino (v^) is disfavored by the kinematics and the G constraint. 
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4.1 .1 Decay of the Inflaton to the Lightest RH Neutrino 

In this case, we suppose that the Majorana masses of vf are hierarchical, with M^c <C M 2 pc, M 3 pc 
(but with M\y C > T r h). The produced lepton-number asymmetry for a normal hierarchical mass 
spectrum of light neutrinos reads [31] 

3 m^Mxpc 

£L = -^^r 5cs - (4 - n) 

Here \S e g \ < 1, which we treat as a free parameter, represents the magnitude of CP violation; m VT is 
the mass of heaviest light neutrino v T and we take (H u ) = 174 GeV - adopting the large tan (3 regime. 

4.1 .2 Decay of the Inflaton to the Next-to-Lightest RH Neutrino 

In this case, we assume M x p a <C M 2 p c <C M 3 p c and impose the conditions T r h < M 2 p c < mi/ 2 and 
M\yc > T r h- The resulting lepton asymmetry is [19,32]: 

^ = 3 M 29 c (m|p - m 2 2D ) 2 s|c| sin 26 
8vr M 3 p c (H u ) 2 (m 2 D sl + m| D c|) ' 

where m^D are the Dirac masses of v j - in a basis where they are diagonal and positive - c# = 
cos "d, s$ = sin with -d and 5 being the rotation angle and phase which diagonalize the Majorana- 
mass matrix, Mpc, of Df. The values of the various parameters are estimated at the leptogenesis scale. 
Note that since Mipc > T r h, el calculated by Eq. (4.12) is not washed out due to v\ inverse decays 
and AL = 1 scatterings - the case with M x pc < T r h is treated in Ref. [56]. Also Eq. (4.12) holds 
provided that M 2 pc <C M 3 pc and the decay width of uf is much smaller than (M 2 ~ c — M^ c )/M 2 pc. 
Both conditions are well satisfied here - see Sec. 4.3.2. 

Since recent results [36] - see, also, Ref. [57] - show that the mixing angle #13 can be taken (at 
95% c.l.) equal to zero and earlier analysis [58] of the CHOOZ experiment [59] suggests that the solar 
and atmospheric neutrino oscillations decouple, we are allowed to concentrate on the two heaviest 
families ignoring the first one. This assumption enables us to connect this leptogenesis scenario with 
the oscillation parameters of the — v T system. The light-neutrino mass matrix, m u , is given by the 
seesaw formula: 

m u = -mbM^rriB, (4.13) 

where ttt-d is the Dirac mass matrix of the z^. The determinant and the trace invariance of m\m v imply 
two constraints on the parameters involved, i.e. 

TO 2D m 3D = m UT m u ^M 2 p c M 3 pc and yl D s| + B D s$ + C D = (4.14) 

with the coefficient of the latter equation being 

= (m| D - ™1d) ((™2D - ™3d) M 2^ + (™2D - m 3D) M 3u- 

+ 2(ml D +m 2 3D )M 2 pcM 39 ccos26j (4.15a) 
B D = 2(m 2 2B - ml D )(ml D M 2 Vc - m 2 m M 2 p c 

2 , _2 



- (™-2D + m 3D )M 2 pcM 3 pc cos 25 j (4.15b) 

C D = {mi B Ml~ c +mj B Ml Pc )-(ml T +ml^ M^M^. (4.15c) 

Assuming that the Dirac mixing angle (i.e., the mixing angle in the absence of RH neutrino Majorana 
masses) is negligible, we can identify [32] the rotation angle which diagonalizes m v with the physical 
mixing angle in the — v T leptonic sector, Q 23 , which is constrained by the present neutrino data - 
see below. 
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4.2 Imposed Constraints 

The parameters of our model can be further restricted if, in addition to the inflationary requirements 
mentioned in Sec. 3.2 and the restriction A > 0.001 which allows us to ignore effects of instant 
preheating - see Appendix B - we impose extra constraints arising from the post-inflationary scenario. 
These are the following: 

4.2.1 To ensure that the inflaton decay to RH neutrinos is kinematically allowed we have to 
impose the constraint: 

mi > 2M K c => M IP c < Xmp/AVScK ~ 1.5 • 10 13 GeV for A > 1.3 ■ 10~ 3 , (4.16) 

where we make use of Eq. (4.2). More specifically, we require mi > 2M\ pc \m\ > 2M 2 ^ C ] for the 
scenario described in Sec. 4.1.1 [Sec. 4.1.2]. 

4.2.2 In agreement with our assumption about hierarchical light neutrino masses for both 
inflaton-decay scenaria, the solar and atmospheric neutrino mass squared differences Attt-q and Am^ 
can be identified with the squared masses of and v T , m% and to 2 t , respectively. Taking the central 
values of the former quantities [36], we set: 

(a) = ^Am 2 Q = 0.0087 eV and (b) m„ T = ^Am| = 0.05 eV. (4.17) 

We multiply the values above by 1.12 in order to roughly approximate renormalization group effects 
for the evolution of these masses from the electroweak up to the leptogenesis scale - see Fig. 4 of 
Ref. [60]. The resulting m VT [T r ^) is low enough to ensure that the lepton asymmetry is not erased by 
lepton number violating 2 — > 2 scatterings [61] at all temperatures between T r h and 100 GeV. Also 
sin 2 623 is to be consistent with the following 95% c.l. allowed range [36]: 

0.41 < sin 2 fl 23 < 0.61. (4.18) 

4.2.3 Due to the presence of S77(4)c symmetry in Gps, m,3D coincides with the value of the 
top quark mass, m t , at the leptogenesis scale - see also Eq. (2.6) - i.e. 

m 3D (T rh ) =m t (T rh ). (4.19) 

Working in the context of the MSSM with tan j3 = 50, and solving the relevant renormalization group 
equations, we find mt(T r h) ~ (120 — 126) GeV for the values of T r h encountered in our set-up. 

4.2.4 The implementation of BAU via non-thermal leptogenesis dictates [14] at 95% c.l. 

Y B = (8.74 ± 0.42) • 10~ n =► 8.32 < Yb/10 -11 < 9.16. (4.20) 

4.2.5 In order to avoid spoiling the success of the BBN, an upper bound on Y~ is to be 
imposed depending on the G mass, mg, and the dominant G decay mode. For the conservative case 
where G decays with a tiny hadronic branching ratio, we have [35] 
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for mg ~ < 



0.45 TeV 
0.69 TeV 

(4.21) 

10.6 TeV 
13.5 TeV. 



The bound above can be somehow relaxed in the case of a stable G. As we see below, this bound 
is achievable within our model model only for > 10 TeV. As gets larger than this bound 
leads to the necessity of a rather fine tuned SUSY breaking mechanism such as split SUSY [63], or 
anomaly mediated SUSY breaking [64] where the superpartners of SM particles have masses lower 
than rriQ. Using Eq. (4.10) the bounds on Y~ can be translated into bounds on T r h- Specifically we 
take T rh ~ (0.53 - 5.3) • 10 8 GeV [T rh ~ (0.53 - 5.3) • 10 9 GeV] for Y d ~ (0.1 - 1) • 10~ 13 

G 



[Y d ~ (0.1 - 1) • 10 -12 ]. 
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Figure 2: Allowed (shaded) regions in the X — Mipc plane, for Xh = Xfj = 0.5, ks = 1 and 
2/33 = 0.5 when the inHaton can decay into v\ 's. The conventions adopted for the various lines and 
shaded or hatched regions are also shown. 

4.3 Numerical Results 

Considering the constraints above in conjunction with those quoted in Sec. 3.2, we can delineate the 
overall allowed parameter space of our model. Recall that we set \h = = 0.5, ks = 1 and g = 0.7 
with Mps given by Eq. (3.27). The upper [lower] bound of the used A's comes from Eq. (3.29) [the 
conclusions of Appendix B]. Also we fix 2/33 = 0.5 throughout our investigation. As can be deduced 
by Eqs. (4.8) and (4.9), for lower [larger] 1/33 's satisfying Eq. (4.20) requires slightly [lower] larger 
Mjpc's. As mentioned in Sec. 4.1, we adopt two leptogenesis scenaria depending on whether Sh can 
decay to v\ or As we see below, in both cases, two disconnected allowed domains arise according 
to which of the two contributions in Eq. (4.3) dominates. The critical point (A c , c^ c ) is extracted from: 

1 - I2c nc ml s = cue = l/12m| s or A c ~ 10" 4 /25.2m| s ~ 0.06 (4.22) 

where we make use of Eq. (3.23) in the last step. From Eqs. (4.7), (4.8) and (4.9) one can deduce 
that for A < A c , T r h remains almost constant; Y\ v c/Y\ decreases as c-ji increases and so the Mjpc's, 
which satisfy Eq. (4.20), increase. On the contrary, for A > A c , T^c /Y\ is independent of c-ji but T r h 
increases with en and so fulfilling Eq. (4.20) Mj^c 's decrease. 

In the following, we present the results of our investigation in each case, separately. 

4.3.1. Decay of the inflaton into v\. In this case our cosmological setting depends on the fol- 
lowing parameters: 

A, cn, M lP c and S eS . 

Given our ignorance about <5 e ff in Eq. (4.11), we take <5 e ff = 1, allowing us to obtain via Eq. (4.11) the 
maximal [62] possible \el\- This choice in conjuction with the imposition of the lower bound on Yb 
in Eq. (4.20) provides the most conservative restriction on our parameters. 

As we explain in Sec. 3.3, M\pc does not crucially influence the inflationary observables. On 
the contrary, M^c plays a key -role in simultaneously satisfying Eqs. (4.16), (4.20) and (4.21) - see 
Eqs. (4.8) and (4.10). For this reason we display in Fig. 2 the allowed regions by all the imposed 
constraints in the A — Mipc plane. In the horizontally lined regions Eq. (3.28a) holds, whereas in the 
vertically hatched region we get h* > mp. The restrictions on the parameters arising from the post- 
inflationary era are depicted by solid, short dotted and dot-dashed lines. Namely, the solid [dot-dashed] 
lines correspond to the lower [most conservative upper] bound on Yb [Yx] in Eq. (4.20) [Eq. (4.21)]. 
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Figure 3: Allowed values of M 2 v versus X, for the input quantities listed in the table above, Xh = 
Xfj = 0.5, ks = 1 and t/33 = 0.5 when the inhaton can decay into 9§ 's. The conventions adopted for 
the types and the color of the various lines are also shown. 



Since we use \5 e g\ = 1, it is clear from Eqs. (4.8) and (4.11) that values of Mypc above the solid line are 
compatible with the current data in Eq. (4.20) for conveniently adjusing \5 c s\ < 1. However the values 
of Mi$c can be restricted by the bounds of Eq. (4.21). Specifically we obtain Yq ~ (0.1 — 1) • 10~ 12 
\Yq ~ (0.5 — 1) • 10~ 13 ] in the [gray] dark gray area. On the other hand, the kinematical condition 
depicted by a short dotted line - see Eq. (4.16) - puts the upper bound on M^c in a the upper right 
corners of the allowed region. As anticipated above, this region has two disconnected branches. In the 
left [right] branch, the first [second] term in Eq. (4.3) dominates. All in all we obtain 

0.39 < M ipc /10 n GeV < 154 for 0.001 < A < 0.062, (4.23a) 
154 > M ipc /10 n GeV > 4.32 for 0.062 < A < 0.25. (4.23b) 

The overall minimal [maximal] M\yc can be found in the left lower [upper right] corner of the allowed 
region. Obviously, the maximal allowed M^c ~ m\/2 is obtained for 0.02 < A < 0.062 and 
0.062 < A < 0.25. This region gives also a lower bound on |<5 c fr|> l^eflfl ^ 1-6 • 10~ 4 . 

Trying to compare, finally, the resulting allowed parameter space here in the A — M^c plane with 
the one allowed within the models of HI [15] - see Fig. 2 of Ref. [65], where the coupling constant 
k corresponds to A - we remark that in our case (i) higher A's and M^c's are allowed and (ii) there is 
an additional minor slice of allowed parameters exclusively due to the SUGRA induced decay channel 
found in Eq. (4.3). 

4.3.2. Decay of the inflaton into v 2 - In this case, our cosmological setting depends on the 
following input parameters: 

X, en, M 2 y C , M 39 c and S. 
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In our numerical program, for every A and c-r consistent with the inflationary requirements of Sec. 3.2, 
we can resolve Eq. (4. 14) w.r.t ?rt2D and if we use M 2 pc and S as input parameters - recall that msB 
is determined by Eq. (4.19). Diagonalizing m v and employing Eqs. (4.8) and (4.12) to estimate Yb, 
we can restrict the parameters through Eqs. (4.18) and (4.20). In order to compare the value of sin 2 #23 
extracted at the leptogenesis scale, T r h, with the low energy experimental result of Eq. (4.18), we solve 
the relevant renormalization group equations following Ref. [60]. We remark that sin 2 #23 increases 
by almost 8% due to these renormalization effects. 

In Fig. 3 we display the allowed values of M 2 „c versus A for various M 3 p c and <5's indicated on 
the left upper corner of the graph. The obtained allowed ranges of several other quantities involved are 
arranged in the Table below the graph. Along the displayed curves the central value of Yb in Eq. (4.20) 
is achieved and we obtain T rh ~ (0.7 - 2.9) • 10 9 GeV which is translated as Y~ ~ (1.3 - 5.5) • 10~ 13 
via Eq. (4.10). From our results we observe that increasing M 3 p c entails an increase of M 2 pc too 
whereas —S approaches tt/4 - see Eq. (4.12) - and increases as M 2 pc drops in order Eq. (4.18) to be 
met. The resulting ?rt2D turns out to be in the range (1.5 — 21.9) GeV which is larger than the mass of 
the charm quark at T r h, m c ~ 0.427 GeV. Therefore within our scheme the SU(4)c symmetry does 
not hold in the up sector of the second family. 

As for the case of Sec. 4.3.1, we obtain two separate branches of allowed parameters, h* > mp for 
A < 0.0037, whereas Eq. (3.28a) is satisfied for A < 0.016 - Eq. (3.28b) is valid for the used A's, see 
Sec. 3.3. Finally, it is remarkable that the assumptions on Mipc, T r h < Mipc <C M 2 pc can be fulfilled 
in a wide range, i.e., 7 • 10 s < M W c/GeV < 10 11 . 

5 Conclusions 

In this paper we attempted to embed within a realistic GUT, based on the PS gauge group, one of the 
recently formulated [10] SUSY models of chaotic inflation with non-minimal coupling to gravity. We 
showed that the model not only supports non-MHI driven by the radial component of the Higgs field, 
but it also leads to the spontaneous breaking of the PS gauge group to the SM one with the GUT break- 
ing v.e.v identified with the SUSY GUT scale and without overproduction of monopoles. Moreover, 
within our model, we can resolve the strong CP and the u, problems of the MSSM via a Peccei-Quinn 
symmetry breaking transition. Inflation is followed by a reheating phase, during which the inflaton 
can decay into the lightest or the next-to-lightest RH neutrino allowing, thereby for non-thermal lep- 
togenesis to occur via the subsequent decay of the RH neutrinos. Although other decay channels to 
the SM particles via non-renormalizable interactions are also activated, we showed that, in both cases, 
the production of the required by the observations BAU can be reconciled with the observational con- 
straints on the inflationary observables and the G abundance, provided that the (unstable) G masses 
are greater than 10 TeV. In the first inflaton-decay scenario, we restrict the lightest RH neutrino mass 
to values of the order (10 11 — 10 13 ) GeV whereas in the second scenario, extra restrictions from the 
light neutrino data and the SU (4)c factor of the adopted GUT gauge group can be also met for masses 
of the heaviest [next-to-lightest] RH neutrino of the order of 10 15 GeV [(10 11 - 10 13 ) GeV] . 

Finally, we would like to point out that, although we have restricted our discussion on the PS 
gauge group, non-MHI analyzed in our paper has a much wider applicability. It can be realized within 
other GUTs, which may (as in the case of Gps) or may not lead to the formation of cosmic defects. 
If we adopt another GUT gauge group, the inflationary predictions and the post-inflationary evolution 
are expected to be quite similar to the ones discussed here with possibly different analysis of the 
stability of the inflationary trajectory, since different Higgs superfield representations may be involved 
in implementing gauge symmetry breaking to Gsm _ see > e -g-> R e f- [66] which appeared when this 
work was under completion. 
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Appendix A: Derivation of the Mass Spectrum During non-MHI 

In this Appendix, we describe the derivation of the mass spectrum of the model when the radial com- 
ponent, h, of the fields v c H , 9 C H slowly rolls down Vui, breaking Gps down to Gsm- We explain the 
results summarized in Table 2, working exclusively in the EF. We demonstrate below the origin of the 
masses of the scalars (Sec. A.l), gauge bosons (Sec. A.2) and fermions (Sec. A.3). 

A.1 Masses for the Scalars 

Expanding Vm in Eq. (3.12) to quadratic order in the fluctuations around the trajectory in Eq. (3.4), 
for given h, we obtain 



V m = Vhio + \™§S 2 + 2 
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where x = u,e,d and g and the decomposition of the scalar fields into real and imaginary parts is 
shown in Eq. (3.7). The various mass-squared matrices involved in Eq. (A.l) are found to be 
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u and e and mi , mi, mi and also mi are presented in Table 2. The elements of the 

6+ 9 g s r 



remaining matrices above are found to be 



mi, = mj,xi(3g'f + A^^)/24/^, mi, = -m^ h {\\x{ - 6) + 3^/)/24/ 
m\ = m 2 P xl(\ 2 xl + 24A^/)/24/ 2 , mi = -X 2 mlx\(x\ - 6)/24/ 2 , 



(A. 3a) 
(A.3b) 



whereas mi g is the same as mi but with Xh replaced by A^-. To simplify our formulae below, we 

take Xh — X H . The various masses squared in Eqs. (A.2a) and (A.2b) originate mainly from Vhf in 
Eq. (3.2). Additional contributions from Vhd// 2 in Eq. (2.15) arise for m| and ml,. The orthog- 
onal matrix Uk, which diagonalizes Mk in Eq. (3.9), diagonalizes M| too. We find the following 
eigenstates 

respectively. The second eigenstate corresponds to the Goldstone 

mode absorbed by A 1 - - see Sec. A.2 - via the Higgs mechanism. Upon diagonalization of M 2 X and 
M 2 2 with x = u, e and d, we obtain the eigenvalues 



9± 



(A.4) 



with eigenvalues m| and m| 
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(A.5) 
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which correspond to the following eigenstates respectively: 

xi± 



1 - \ ~ 1 /- ^ n 

— P (xi ± xij and X2=p = — f= (x2 =F with x 
\/2 \/2 



■u a , e and d. 



(A.6) 



Note that ml 



JO 



A rripxi/Af <C m|_ for x = u (3 colors) and a; = e are the masses 
squared of the Goldstone bosons x\ + and a?2- which are absorbed by A^ (3 colors) and A^ - see 
Sec. A.2 - via the Higgs mechanism. The remaining m| ± 's are listed in Table 2. 

A.2 Masses for the Gauge Bosons 

Some of the gauge bosons A^ and A 7 ^ acquire masses from the lagrangian terms - cf . Ref . [39] : 

K a0 {{D^Y (D»H cT f + {D^) a (D^H c f^j = 
\m\ ( (J\A^ - AQ (J*A$* - 4^) + + 2 ELi . (A-7) 

where m± are given in Table 2. The action of on i7 c and is as follows 



cT 



d,H cT + ig ( Eili TqAqH^ + Z' m=1 T™A™H< 



a 5 



and we have defined the following normalized gauge fields: 



A^ = -L (A?+ 2a ± iA^) for a = 1, 2, 3 and A± = ± ± iA^) . 



C ~ V2 K ° ° 
The first term of the RHS of Eq. (A.7) can be written as 



' A 1 ^ 

4" 



m\A^Ai 



(A. 8a) 
(A. 8b) 

(A.9) 
(A. 10) 



where m\ = 5m±/2 - see Table 2 - since 



3/2 
V372 



/372 
1 



and U CK M CK Ul K = diag (5/2, 0) 



with 



^CR = 



3/5 V2/5 
275 









■a^ 


and 


AH 


= U C R 




> 





(A. 11) 



(A. 12) 



Therefore, from Eqs. (A.7) and (A. 10), we can infer that 9 gauge bosons (A^, A^ 1 and A^) become 
massive, absorbing the massless modes 0_, uf + , ei+ and e2_, whereas All, which remains mass- 
less, can be interpreted as the B boson associated with the U{l)y factor of Gsm- 

A. 3 Masses for the Fermions 

The lagrangian kinetic terms of the chiral fermions ip x with x = v, e, u, d, g, v, e, u, d, g, which are 
associated with the superfields Vjj, e° H , u° H , d c H , g c , v c H , e c H , u c H , d c H , g c respectively, are 



= $ v+ ^D$ v+ + $ v _^D$ v - + , £ i $ y ai 1 D li $ v , (A.13) 
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where we have used Eq. (3.9), and the normalized spinors are defined as follows 



l 



with t/v± = —7= ± tyv) and ip y = — % 



(A. 14) 



^ ^v--— — ^ ^ 

with y = e,u,d,g,e,u,d, and g. In Eq. (A. 13) the contraction between two Weyl spinors is sup- 
pressed; the Pauli matrices and the action of on ip a are specified in Ref. [3,42]. 

Having defined the normalized spinors, we can proceed with the derivation of the fermionic mass 
spectrum of our model. The masses of the chiral fermions can be found applying the formula [3,42]: 



m Q(3 



e */2m£ (V HP Sa/3 + \ (^a/jWaPS + K a F fi + Kf}F a ) - T^k) (A.15) 

with = K^daKpy, Wupsa/3 = Wnps^c,^, K a = K^ a and F a as defined below Eq. (3.2). 
Upon diagonalization of the relevant mass matrix, we obtain the eigenvalues , m^_ d and , 
listed in Table 2, corresponding to the following eigenstates 

^s»± = j= ($s ± Vw) , ^ ± = (% ± and = -L ± ^) . (A.16) 



We remark that Whps in Eq. (2.9) does not give rise to mass terms for fermions in the sectors u c H — u c H 
and e c H — e c H . However fermion masses also arise from the lagrangian terms 

-iV2gK a? ( ZlU Ac (fn* (T a H cJ f - Ti (T«*H*f) + 

YL=i ^ (Thc {T m H cJ f - (T™H c f) + h.c.j , (A.17) 

where A^ [A R ], is the gaugino corresponding to the generator Tq [T^] and ipgc, i[>h c represent the 
chiral fermions belonging to the superfields H c , H c respectively. Concentrating on T<i 5 , T^, we obtain 



-iV2g (x£ [i>„ 

igh jjp - ip u ( fS 
2/ V2 I V2 



M K 



rpl5* JJC 

C 

T 15 H cl 



+ A 



15 



M K 



_ T 3* H c 
T R HCl 



+ h.C. 



A c 5 + A| +h.c. 



im±ip u -X ± + h.c. with 



A x 
AH 



Ac 5 
Ar; 



Therefore, we obtain a Dirac mass term between the chiral fermion and the gaugino — iX ± , 
whereas a Dirac spinor composed by the combination of ip v+ and — iA" remains massless and can 
be interpreted as the Goldstino which signals the (spontaneous) SUSY breaking along the direction of 
Eq. (3.4). 

Similarly, focusing on the directions J* +2a with a = 1, 2, 3 and and Tj|, we obtain the mass 
terms 



'2^7 



£ (<A C + - r a X^ + h.c.) + (VU R - ^e-A+ + h.c.) , (A.18) 



va=l 



where we have defined the following combinations of gauginos 

A a c ± = ± ^ ± i^) for a = 1, 2, 3 and A± = ± (fa ± zA|) 



(A. 19) 



V2 VV c ' ' ' ^2 

in agreement with the definition of the corresponding gauge bosons in Eq. (A.9). Therefore, the chiral 
fermions ^ and ^ [ip e and combine with A^ [A R ] to form two Dirac (or four Weyl) fermions 
with mass m± as one deduces from Eq. (A.18). This completes the derivation of the spectrum of the 
model along the inflationary trajectory of Eq. (3.4). 
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Appendix B: Inflaton Oscillations After non-MHI 

In this Appendix we discuss various (p)reheating mechanisms [51] which could become competitive 
with the perturbative decay of the inflaton to lighter degrees of freedom, as analyzed in Sec. 4. Indeed, 
in certain regions of the parameter space, the process of reheating in this theory can be quite com- 
plex. After the end of non-MHI, the inflaton develops a tachyonic mass, crosses an inflection point 
and enters into a phase of damped oscillations. As pointed out in Ref. [52], where a similar potential 
is investigated, the particle production due to tachyonic preheating is not significant because the pas- 
sage of the inflaton through this region is very short. During the subsequent oscillations, perturbative 
production of superheavy bosons - i.e., bosons with masses at the SUSY vacuum proportional to Mps 
- is not possible, since these particles are heavier than the inflaton at the global minima of its potential 
as shown in Table 3. Therefore effects of narrow parametric resonance [51] are also absent. However, 
if the initial amplitude of the inflaton oscillations is large enough, it may pass through zero, h = 0, 
where these bosons are effectively light and can be produced through (non-perturbative) instant pre- 
heating [53], as we discuss in Sec. B.2 below. We first study the dynamics of the inflaton' s oscillations 
in Sec. B.l. 

B.1 Dynamics of the Inflaton Oscillations 
The cosmological evolution of h (h) in the EF is governed by the equation of motion: 

h + 3Hh + V mQ £ = (B.l) 

where H is the Hubble parameter in the EF and Vhio is given in Eq. (3.5) - recall that the dot denotes 
derivation w.r.t. the cosmic time t and h = Re(z^ + v c H )/2 along the direction of Eq. (3.4). In the LHS 

of Eq. (B.l), we neglect the damping term T\h which is important only at the stage of rapid oscillations 
of h near one of the minima of Vhio [51]- Note that contrary to the case of the potential analyzed in 
Ref. [52], the minima of Vhio lie at = 2Mps S> and Vhio nas a maximum at h = with 

Vmo(h = 0) = Vb = A 2 M| S (B.2) 

which can not be ignored. Due to these features, the quadratic approximation to Vhio is not accurate 
enough for the description of the h post-inflationary evolution. 

The solution of Eq. (B.l) can be facilitated if we use as independent variable the number of e- 
foldings N defined by 

N = In (r/R^ => N = H and H = H'H . (B.3) 

Here the prime denotes derivation w.r.t. N, R(t) is the EF scale factor and R\ is its value at the com- 
mencement of non-MHI, which turns out to be numerically irrelevant. Converting the time derivatives 
to derivatives w.r.t. N, Eq. (B.l) is equivalent to the following system of two first order equations 

F h = JHR 3 ti and JHF' h = -Vhio,/^ 3 with F h = hR 3 . (B.4) 
This system can be solved numerically by taking 

H = -j^— (Fl/2R & + Vhio) ^ (B.5) 

and imposing the initial conditions (at N = 0) h(0) = (0.5 — 2.5)mp and h'(0) = 0. We checked that 
our results are pretty stable against variation of h(0). 
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Figure 4: The evolution of the quantities h/mp (upper plots) and \m±/rn±\ (lower plots) as func- 
tions ofN~ N{ for A = 0.0037 and c n = 81 (left panel) or A = 0.01 and c n = 235 (right panel). 



During non-MHI we have H ~ -ffmo and the results of Eqs. (3.20) and (3.23) are well verified. 
Soon after the end of non-MHI, we obtain H ~ Huwe~ 3 ( N ~ Nl ^ 2 - with N t being the value of N at 
the end of non-MHI - and h enters into an oscillatory phase with initial amplitude equal to /if given by 
Eq. (3.19). Since the value of Vhio at the end of inflation, Vhio is larger than the value of Vhio at its 
local maximum h = 0, Vo, we expect that h crosses zero at least once during its evolution. However, 
as can be deduced from Eqs. (3.19) and (3.23), lowering A increases /if but decreases Vo. Therefore 
the passage of h through zero is facilitated. 

The intuitive results above can be established and refined through the numerical solution of Eq. (B .4), 
during the h oscillations, depicted in the upper plots of Fig. 4. Namely in left [right] plot we present 
the evolution of h as a function of N — TVf for A = 0.0037 and c-ji = 81 [A = 0.01 and c-ji = 235]. In 
both cases, we see that h, decreasing slowly from its value /if = 0.152 [/if = 0.056] for A = 0.0037 
[A = 0.01], passes from the minimum of Vhio at /i = 2Mps and then climbs up the hill of Vhio at 
h = 0, falls towards the other minimum of Vhio at h = — 2Mps until it reaches a maximal value and 
oscillates backwards. This path is followed some times until h falls finally into one of the minima of 
Vhio at h = — 2Mps [h = 2Mps] for A = 0.0037 [A = 0.01] - performing damped oscillations about 
it. In other words, h oscillates initially around the local maximum of Vhio and then about one of the 
two SUSY vacua. The number of passages though zero increases as A decreases - it is equal to 4 [12] 
for A = 0.01 [A = 0.0037]. Solving repetitively Eq. (B.4) we notice that h ceases to cross h = for 
A > 0.088. 

B.2 Instant Preheating 

Whenever h crosses zero particle production may occur via instant preheating [53]. This mechanism is 
activated when the /i-dependent effective masses, m e s, of the produced particles violate the adiabatic- 
ity criterion, according to which 



l^cfr/^cfrl = Hm' eS /m 2 eS 



< 1. 



(B.6) 



Here, m c g represents collectively the masses of superheavy bosons with masses proportional to <?Mps, 
A#Mps or XgMps - see Table 2. We focus on the production of these bosons since these can sub- 
sequently decay efficiently to the light SM particles altering drastically the picture of the usual pertur- 
bative reheating. On the contrary, the scalars of the S — v c H — v c H sector with masses proportional to 



AMps have suppressed decay modes to the RH neutrinos only. Taking as an example m e ff = m± we 
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plot in the lower plots of Fig. 4 the evolution of |m±/m±| as a function of N — Nf for A = 0.0037 
and c-r. = 81 [A = 0.01 and en = 235] - see left [right] plot. We observe that Eq. (B.6) is violated 
more frequently as A drops since the passages of h through zero become also more frequent. From 
the results of our numerical treatment we find that Eq. (B.6) holds during the whole post-inflationary 
evolution of h for A > 0.045 whereas it fails more than 40 times for A < 0.001. 

The produced this way superheavy bosons acquire a large mass while the inflaton increases towards 
its maximum amplitude and start to decay into all lighter particles within almost half oscillation of the 
inflaton, rapidly depleting their occupation numbers. As shown in Ref. [52] an efficient transfer of 
energy from h to the superheavy bosons requires a rather large (let say 50 — 70) number of passages 
of h through zero. Meanwhile effect of backreaction of the produced particles on the h condensate 
may become significant and a more involved numerical study of the process is imperative. Trying to 
deliberate our leptogenesis scenario from such a complicate situation, we impose the indicative lower 
bound A > 0.001 - which can be translated as a bound c-r > 21 via Eq. (3.23) - above which our 
estimations in Sec. 4 are more or less independent of the preheating effects. We hope to return to the 
analysis of A < 0.001 region in the future. 
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